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We investigate the antiadiabatic limit of an antiferromagnetic 5=1/2 Heisenberg chain coupled to 
Einstein phonons via a bond coupling. The flow equation method is used to decouple the spin and the 
phonon part of the Hamiltonian. In the effective spin model longer range spin-spin interactions are 
generated. The effective spin chain is frustrated. The resulting temperature dependent couplings 
are used to determine the magnetic susceptibility and to determine the phase transition from a 
gapless state to a dimerized gapped phase. The susceptibilities and the phase diagram obtained via 
the effective couplings are compared with independently calculated quantum Monte Carlo results. 

PACS numbers: 63.20.Ls, 75.10.Jm, 75.50.Ee, 63.20.Kr 
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I. INTRODUCTION 

We study the magnetic susceptibility and the zero- 
temperature phase diagram of the antiferromagnetic S = 
1/2 Heisenberg chain coupled to dispersionless Einstein 
phonons. Due to the coupling of the lattice degrees of 
freedom to quasi one-dimensional magnetic degrees of 
freedom, this model exhibits a zero-temperature phase 
transition between a gapless phase and a massive phase 
showing dimerization. This quantum phase transition 
in the one-dimensional model is regarded as the equiv- 
alent to the spin-Peierls transition!!! occurring at finite 
temperature in the three-dimensional model. Models 
of this type wete investigated earlier by Pyttetl and 
Cross and FisherB. The discovery of the first inorganic 
spin-Peierls substance CuGe0 3 aE3 renewed the interest 
in models with spin-phonon (SP) coupling. n Various 
methods (e. g. densitjs-, matrix renormalizationQ, xutoff 
renormalization g rn |jry, linked cluster, pypansinrfi ex- 
act diagonalizatii 



flow equationsE3E3, and quan- 
tum Monte CarloE^t-JtL3) were used to investigate ground 
state properties, especially the phase diagram, the mag- 
netic excitation spectrum and thermodynamic properties 
of the spin-phonon chain. 

In this paper we compare results of two independent 
methods applied to the spin chain Hamiltonian with 
a local bond coupling to phononic degrees of freedom: 
flow equations (continuous unitary transformations) and 
quantum Monte Carlo. With the flow equation method 
we map the initial magneto-elastic problem onto an ef- 
fective magnetic problem with temperature dependent 
effective spin-spin couplings. These are used to calculate 
the magnetic susceptibility (via high temperature series 
expansions or exact complete diagonalization) and to de- 
termine the T — phase separation line. Independently, 
the phase diagram and x(T) are determined by means of 
quantum Monte Carlo (QMC). 

In Sect. II we introduce the model Hamiltonian for 
the spin-phonon chain and present the flow equation ap- 
proach used. Next, the quantum Monte Carlo approach 



is introduced in Sect. [II. Finally, we show our results 
for the T = phase diagram in Sect. [V A and for the 



IV B. A short summary 



magnetic susceptibility in Sect, 
is given m Sect. 



II. DEFINITION OF THE MODEL AND THE 
FLOW EQUATIONS 

The model under study reads 

H{£) =H S + Hsb + H B (la) 
Hs = ^(JiWS 4 S. i+1 + J 2 (£)S l S l+2 ) (lb) 

i 

H B =uJ2 b i bi ( lc ) 

i 

H SB = J2(Mt)b\+4(e)h) ■ (Id) 

i 

Here Si stands for the 5 = 1/2 spin operator on site i and 
bi destroys (creates) a phonon on the bond between 
the sites i and i + 1. The variable t is used to parame- 
terize the flow equation transformation and is explained 
in detail in the following paragraph. Tia, is the Hamil- 
tonian of a frustrated spin chain or J\-Ji model. If the 
frustration parameter a = J2/J1 exceeds a critical value 
of a — a c = 0.241167(5), the model undergoes a quan- 
tum pbasfi— transition from a gapless state to a gapped 
phaseN'tliy. We start with the Hamilton operator (Q) 
consisting of pure spin (Ws) and phonon (Wr) parts and 
a spin-phonon coupling term Hsb- Flow equations are 
used to decouple the spin-phonon system. 

The flow . equ ation approach was introduced by Weg- 
ner in 1924E3E3. This method is similar to Frohlich's 
approachEiJ, which maps a Hamilton operator onto an ef- 
fective Hamiltonian in one step by applying a unitary 
transformation. Wegner's approach constitutes a modi- 
fication of Frohlich's single step transformation since an 
infinite number of infinitesimal unitary transformations 



2 



is used. The idea behind this scheme is that the contin- 
ued adjustment of the infinitesimal transformation to the 
Hamiltonian yields a smoother effective interaction than 
a rotation in one step. 

The continuous unitary transformation is parameter- 
ized by a flow parameter I £ [0, oo]. So H(£ = 0) stands 
for the original Hamiltonian. The effective Hamiltonian 
7i{£ — oo) is simpler in the way that the direct spin- 
phonon coupling has been rotated away. The infinitesi- 
mal generator rj(£) defines the unitary transformation via 



dH(£) 
d£ 



W),H{1)] 



A choice of rj{£) proposed by Wegner is 
V (l) = [H d (l),H(£)} , 



(2) 



(3) 



where Tid is the part of the Hamiltonian which is taken 
as the diagonal part. By choosing Hd(£) = 'Hs{£)+T~Lb{£) 
and the "off-diagonal" interaction part as Ti d{£) = 
Wsb(^), our generator reads n(£) = [Wd,W d] = [Hs + 
7~Lb , HsbIh i AJternative choices of rj were used in the 
literatureEjSEa in order to obtain flow equations with 
specific properties. For the present purposes, the canon- 
ical choice (||) of r\ is sufficient. 
As initial conditions we choose 



J x (0) = J and J 2 (0) = J 2 ° = 0. 



(4) 



So the untransformed pure spin part is the antiferromag- 
netic S — 1/2 Heisenberg Hamiltonian with a nearest 
neighbor coupling. In the limit £ — > oo longer range 
spin-spin interactions arise and the next-nearest neighbor 
coupling constant J 2 m the effective spin model becomes 
finite. 

Ti.SB describes the coupling between spin and phonon 
system. The coupling operator Ai(0) typically consists 
of nearest neighbor spin products since in realistic mag- 
netic materials the spin-phonon coupling mainly influ- 
ences neighboring sites. By choosing different coupling 
operators various mechanisms of how the lattice distor- 
tions influence the exchange integral J can be investi- 
gated. If the exchange coupling depends on the position 
of a spin between the neighboring ones the appropriate 
choice for Ai(0) is the difference coupling 



Af lS — g(SiSi+\ — S;S^i) . 



(5) 



For CuGeC>3 this corresponds to changes of two neigh- 
boring Cu-O-Cu binding angles induced by shifts of the 
copper ion in the middle, i. e. one angle is enlarged at 
the expense of the other one. This coupling type was in- 
vestigated in the framework of a flow equation approach 
in Refs. ll,lj. In Rcf. |ll| an analysis is presented in 
leading order in g/io and in the two leading orders in 
J/uj. Though only systematic to order g 2 /uj 2 and g 2 J/uj 3 
a good qualitative picture for the antiadiabatic limit of 
the spin-phonon chain could be achieved. A more sub- 
tle analysis is given in Ref. [l2| The price paid in Ref. 



|l2| for high precision results for small J/u> and g/ui is a 
breakdown of the formalism for large values of J/uj. 

In this paper we study the bond coupling 



A(0)=5SiSi +1 



(G) 



which is frequently usedBEHHi3ll3'E£l, especially in connec- 
tion with CuGeC>3. For this type of coupling single har- 
monic degrees of freedom directly modify the magnetic 
interaction. This is due to. akin, group effects engendered 
by the germanium atomsoLjEj. In the present paper we 
restrict ourselves to an analysis in the leading orders as 
in Ref. 

To study the Hamiltonian (Q) with the coupling opera- 
tor (Q) via the flow equation method, let us first observe 
that according to Ref. Ai should be normal-ordered, 
i.e. Ai(0) — > Ai(0) — {Ai(0)}. In this way the quan- 
tum phonons only couple to the comparably small spin 
fluctuations and not to the bare operator S^Si+i. This 
is an extremely important point. Only the coupling to 
the small fluctuations justifies the expansion approach. 
This is not guaranteed by the specific choice (|6j). This 
feature of the coupling type (|^) is in contrast to the dif- 
ference coupling Af lS which does not contribute to the 
undistorted phase due to the translational invarianceE-a, 
i.e. (Af lS (0)) vanishes. The necessary normal-ordered 
coupling operator for the bond coupling reads 



A(0) = g(SiS i+1 - C) 



(7) 



where C stands for the (temperature dependent) expecta- 
tion value of the nearest neighbor spin product (SiSi + i). 
The value of C is calculated using Pade approximants for 
the series expansions given in Refs. p4p9| . As we deter- 
mine C in the uniform phase it can be interpreted as an 
average value over the chain For the antiferromagnetic 
chain the expectation value C for T = lies between 
1/4 - ln(2) w -0.443 (for a = 0) and 3/8 = -0.375 (for 
the fully dimerized chain). 

To achieve normal-ordering we can rewrite the Hamil- 
tonian (Q) by shifting the phonons operators via bi — > 
b\ — —. Thus, the Hamiltonian with the normal-ordered 
coupling is 



■ CO 



J = J — and Eq = — 



LO 



LO 



(8a) 



(8b) 



We now follow the outline of Ref. |Tl|. Applying the uni- 
tary transformation generated by 77, we find an effective 
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Hamilton operator without spin-phonon coupling terms 



H eS = j^SiS i+1 + AHx + AH Y 

i 

i i 

AHx = --J2 A t-Ai 



uj 

.,2 



(9a) 



(9b) 3 



-, w £ K 1 + 4C ) S * S *+i - - 3 /8] 

i 

^2COth(^)^[^,[W S ,A]] 

l^coth^^^-J)^ 



(9c) 



As no frustration is present in the original model (i. c. 
J2 = 0) we obtain for the effective nearest and next- 
nearest neighbor couplings 



Ji = J + t 



9_ 

2uj 



^coth(^) 
2uj 2 \2TJ 



Jg 2 



(10) 



and j 2 = g_ coth (^) 

while longer range two-spin couplings and products with 
four different spins are omitted. In this way the dressing 
of the spins with phonons induces a frustration a e s = 
J2/J1 > 0. The omission of the longer range and the 
four-spin couplings contributes an additional approxima- 
tion. It is justified by the fact that these terms contribute 
neither in a Neel state nor in a dimer state. We note that 
the couplings are temperature dependent via the coth 
terms and via the T-dependence of the expectation value 
C = (SjSj+i). Though only systematic to order g 2 /ui 2 
and g 2 J/uj 3 , higher order effects enter due to the appear- 
ance of J = J — 2 ~ in Eqs. (fio|). Since C is calculated 
self-consistently, this is in principle a correction of infi- 
nite order. In this way the phonon-shift bi — > bi — — , 
which is necessary to apply the flow equation approach, 
induces infinite order contributions. ...... 

It was pointed out in recent articlesOoEJ that the 
spin-Pcicrls transition in the antiadiabatic limit is not 
accompanied by a softening of the Peierls-active phonon 
modes as it is the case for organic spin-Peierls com- 
pounds. This phonon hardening was observed for the dif- 
ference coupling using a similar flow equation approach in 
Ref. For the bond coupling investigated here a sim- 
ple formula for the effective phonon frequency u> c s can 
be derived by evaluating the btb^ term of dAH/d£ (cf. 
formula (9) of Ref. O). Hence we find 



9 2 J 



HSiSi+i) + (SA+a)) (11) 



where non-local terms are neglected. As for zero tem- 
perature (SiSi + i) < 0, (SiSi+2} > 0, and thus J = 




FIG. 1: Temperature dependence of the effective phonon fre- 
quency. The plot shows the relative change in u for cu/J — 4 
and g/J = 1.5. 



J — 2g 2 (SjS i+ i)/cj > J the effective phonon frequency 
is enlarged. With increasing temperature the correla- 
tions are destroyed and accordingly w e ff tends to u. The 
generic temperature dependence of the relative frequency 
increase is depicted in Fig. [l] for uj/ J = 4 and g/J = 1.5. 

After these remarks regarding the phonon part of the 
effective Hamiltonian, we now turn back to the physics 
in the spin sector. By solving a = J2/J1 = a c for g/u> 
for T = we obtain an explicit formula for the phase 
separation line 



7c - \/{lc- a c ) 2 - 16Ca c 7c - Qc 
4C(l+a c ) 



(12) 



with the abbreviation 7 C = (l+a c ) J/uj. The two limiting 
cases are then given by 



lim — 

J/w^o UJ 

lim — 

J j uj — >oo UJ 



2C(l + a c ) 

0.6234 . 



0.4682 



2a c 



1 + a c 



(13a) 
(13b) 



The finitcness of the critical spin-phonon coupling g c /to 
for J /to — > is a striking difference between the local 
coupling mechanism considered here and the coupling to 
the order parameter used e. g. in Ref. |lj. A qualitative 
understanding can be achieved by considering the shifted 
and rescaled Hamiltonian 



H 



J 



2 (5) : 2 c J2 SA + i+ J2^ 



(14) 



If g/10 is finite the effective nearest neighbor coupling 
constant remains finite even in the limit J /uj — > 0, as the 
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FIG. 2: Dependence of the effective nearest and next-nearest 
neighbor spin-spin couplings from the spin-phonon coupling 
g for fixed values of u> and zero temperature. 



expectation value C is typically rs —0.4. Thus, there is 
no reason why the critical spin-phonon coupling should 
vanish for J/ui approaching zero. This illustrates an im- 
portant difference between the difference coupling (^) and 
the bond coupling (H). 

Fig. H shows the change of the effective spin-spin cou- 
plings with increasing spin-phonon coupling. Since our 
method corresponds to the leading orders of an expan- 
sion in J /to and g/to, the approach breaks down if the 
ratio J/uj becomes too large. This can be seen in the 
unsystematic downturn of the Ji(g) curve for lu = 2J. 

To obtain the magnetic susceptibilities for the spin- 
phonon chain we follow various routes. Using the effec- 
tive couplings ( JTo) ) we perform an exact complete diag- 
onalization for the frustrated spin|-chain or use alterna- 
tively high temperature expansionsEH In order to extra u- 
olate the results we use Dlog-Pade approximantsEjEi 
If the spin-phonon coupling g/u> is small, the effective 
frustration is negligible. This justifies the third ap- 
proach to simply rescale the,,pxact results for the un- 
frustrated Heisenberg chai rOT— lO by the renormalizcd 
nearest- neighbor coupling J\ . The results for these three 
methods agree perfectly well with the direct QMC data, 
which are not based on the derived effective couplings. 



Thus our results are not hampered by too low cutoffs in 
the phonon number, which represent otherwise a serious 
difficulty at low phonon frequencies. As has been demon- 
strated in Ref. |l5j the error originating from the finiteness 
of the Trotter number can be easily overcome by extrap- 
olation. Here Trotter numbers of M = 40, 60 and 80 
are used. The calculations are performed with 10 5 -10 6 
spin updates and 20 phonon updates per spin update. 
We considered up to 256 sites so that the finite-size ef- 
fects are controllably small. Thus, apart from statistical 
fluctuations the Monte Carlo results are exact. 

It was appropriate for the quantum Monte (Carlo study 
to use the Hamiltonian in the following form£L3 

n = 2jJ2iSiS i+1 -h[l + ~g(bl + b i )]+Qj2 b l b i ( 15 ) 



where we use the tilde to distinguish clearly between the 
different coupling constants in the different ways to de- 
note the Hamiltonian. The operator ( |l5|) can be mapped 
onto (Q) by a shift in the phonon operators 



bi — > bi 4 

By this mapping one obtains 



Jg 

2u> 



(16) 



n =2j(i + iL)^2 SiSi+i + q b\bi 



(17) 



Accounting for this rescaling of the spin-spin coupling 
constant the susceptibilities computed for H in Eq. ( p"5| ) 
and those for TL in Eq. (|]) can be compared. 

The quantum Monte Carlo procedure is per construc- 
tion a method giving results at finite temperature T. But 
since the algorithm allows to approach very low T one can 
also study properties of the ground state, in particular if 
one considers the model in its gapped phase. This was 
used in the following to extract the phase boundaries of 
the model. 



IV. RESULTS 



A. Phase diagram 



III. QUANTUM MONTE CARLO 

In a completely independent approach we study the 
model (Q) directly by a quantum Monte Carlo method 
based on the loop algorithmLj. The technical details of 
the method are described in Ref. 15|. Here we just point 
out again, that the algorithm al ows to incorporate as 
many phonons as desired. We usually choose 50 per site. 



The calculation via flow equations is based on the idea 
that the transition into the ordered phase does not oc- 
cur due to the softening of a phonon but due to the 
tendency of the effective spin model towards dimeriza- 
tion. Phonon-induced frustration above its critical value 
a Q g > a c drives the dimerization. The phase transition 
line is therefore determined by solving 



lim a 



(s))| 



(18) 



5 



2.5 




number of site 



FIG. 3: Zero temperature phase diagram of the spin-Peierls 
antiferromagnetic chain of spins interacting with phonons. 
For small values of the spin-phonon coupling g/u) the sys- 
tem is gap less (A = 0). For large g/uj the system is dimerized 
and has an energy gap (A > 0). The circles (squares) corre- 
spond to QMC data, where (no) dimerization was found for 
T I J = 0.05 and a chain of N = 256 sites. QMC data can only 
be computed to the right of the dot-dashed line (for details 
see main text). The results from Weifie et. al. are taken from 
Ref. 0. 



for fixed values of J with respect to g (cf. Eq. [lj). In 
this way we find the critical spin-phonon coupling g c in 
dependence of the nearest neighbor spin coupling J as 
depicted in Fig. [| 

Using the quantum Monte Carlo method the decision 
whether a pair (J/u),g/ui) represents a point in the un- 
gapped or the gapped phase is reached by investigating 
the expectation values of the local phonon displacements 
(b\ + bi). At low temperatures (b\ + &j) show clear dimer- 
ization patterns if the model is in the gapped phase and 
they are randomly fluctuating in the critical phase of 
the model (see Fig. ||). Of course, the ergodicity of the 
Monte Carlo procedure complicates the identification of 
the ground state in the gapped phase since for long run 
times the ground state fluctuates always between the two 
possible dimerized patterns. However, by averaging over 
not too many configurations (we generally take 10000) 
one can extract well the dimerization patterns from the 
Monte Carlo results. 

For intermediate values of J/oj the flow equation and 
QMC method predict similar values for the strength of 
the critical coupling g c /uj (cf. Fig. |J). There is a cer- 
tain tendency that the QMC data indicates lower values 
for g c jio than the analytical calculation. We attribute 
this discrepancy partly to the approximate character of 
the analytical calculation and partly to the difficulty to 
determine the phase boundar y by QMC precisely. Note 
that due to the phonon shift ( |l6|) the QMC approach is 
restricted to the region on the right hand side of the dot- 
dashed line in Fig. |3[ The QMC data indicate a smaller 
value for g c /u> in the antiadiabatic limit. In the immedi- 



FIG. 4: Phonon displacement averaged over 10 4 confi gura- 
tions for Q — 10 J and g = 0.8 (g = 0.1) for a system size 
of N = 256 sites in the upper (lower) plot. For g — 0.8 the 
dimerization is clearly visible. 



ate vicinity of the dot-dashed line (i. e. for very large oj) 
it is rather hard to distinguish between a fully dimerized 
pattern and partially dimerized structures. Thus, further 
investigations are desirable for the extreme antiadiabatic 
limit of the phase diagram. 

The phase separation line predicted by Weifie et. alx3 
via fourth order perturbation theory lies well above the 
flow equation curve and above the Monte Carlo results. 
This was also observed in Ref. [l2] where in the, case of the 
difference coupling the leading order resultllJ compared 
much better with the DMRG datacl than the fourth or- 
der resultE3. Since with a difference coupling the critical 
spin-phonon coupling g c /uj vanishes as J/uj goes to zero 
the expansion in these parameters becomes exact. Thus, 
the phase separation lines for the difference coupling co- 
incide for small values of J/lu!E1 But as discussed above, 
in the case of a bond coupling g c /oj is non-zero in the 
antiadiabatic limit (J/oj — > 0). Hence, it is not surpris- 
ing that even in the limit of vanishing J/oj the devia- 
tion does not decrease. Ref . 10 predicts even a value of 

0. 



g c /oJ = a/8/3 • a c /(l + 2a c ) » 0.6587 for J/lj 

In Ref. [l4| a QMC analysis of the spin-phonon chain 
with a bond coupling was presented for the parameter 
J/oj = 4. A critical coupling g c /oj — 0.636 ± 0.042 was 
calculated by investigating the staggered spin suscepti- 
bility for various system sizes. This result is in very good 
agreement with g c /oj — 0.613 predicted by the flow equa- 
tion method (not depicted in the phase diagram Fig. |^). 



B. Magnetic susceptibility 

In this Section we consider the influence of the spin- 
phonon coupling on a thermodynamic quantity like the 
magnetic susceptibility x(T). We show that the magnetic 
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FIG. 5: Magnetic susceptibility for g = 0.5J and to = 4J. 



FIG. 6: Magnetic susceptibility for g = 1.5J and lu — 4J. 



susceptibility can be well described by an effective spin 
model. Furthermore, the question whether it is possible 
to describe thermodynamic quantities of the SP chains 
with a static spin model or via a spin model with tem- 
perature dependent couplings J(T) is investigated. The 
exchange integral J of a quantum spin system is often as- 
sumed to be constant. However, in real materials where 
J depends on the actual positions of the magnetic ions 
in the crystal this approximation need not be valid. For 
the spin-Peierls substance CuGeC>3 this point is of special 
interest. The magnetic susceptibility of CuGeC>3 in the 
temperature range above the spin-Peierls transition, can 
be well fitted by a frustrated Heisenberg model. EZrE§E3 It 
is often objected, however, that the agreement might be 
accidental as no interchain interactions and spin-phonon 
couplings are taken into account. This motivates us - 
although the Hamiltonian introduced in Eq. ([!]) is too 
simple to be a realistic model for CuGeC>3 — to study 
in detail the influence of the spin-phonon coupling on 
the thermodynamic observables. Another substance with 
strong magnetOpelastic couplings, (VO)2P207, was inves- 
tigated recentlyEy by a flow equation approach. The au- 
thors found that even for strong spin-phonon coupling the 
influence of phonons on thermodynamic magnetic quanti- 
ties is rather small. A good description of such quantities 
could already be achieved by a static model. 

The flow equation approach maps the spin-phonon 
chain onto a frustrated spin chain with temperature de- 
pendent spin-spin couplings. The corresponding x(^) is 
compared with the susceptibility calculated with QMC 
for the full Hamiltonian (jl^). This enables us to study 
the influence of the spin-phonon coupling on x{T)- Fur- 
thermore, we can discuss the question whether it is pos- 
sible to retrieve the main physical properties of a spin- 
phonon chain by looking at its effective spin model. 

Figures |^-|| show our results for four different values 
of g and w. All parameters are given according to the 
Hamiltonian (111). As a reference the exact result of the 



unfrustrated Heisenberg chain@S@ is depicted in all fig- 
ures. Fig. [5] illustrates that a rather small spin-phonon 
coupling g = 0.5 J leads only to a lowering of the overall 
height of the susceptibility. The frustration in the effec- 
tive model is very small (J2/J1 ~ 0.008 for T = and 
J2/J1 ~ 0.02 for T/J = 5). This enables us to rescale the 
Heisenberg curve with respect to the enlarged coupling 
constant J\. To ensure that the effective frustration is 
really small enough to be neglected we make an addi- 
tional check via the exact and complete diagonalization 
(ECD) data (computed with the temperature dependent 
couplings J\ and J2) for a frustrated chain with N = 16 
sites. The QMC and the flow equation results agree per- 
fectly down to very low temperatures. 

When g is increased the maximum decreases as can be 
seen in Fig. || for u> = 4 J and g — 1.5 J. Again, the two 
approaches yield the same results down to T/J « 0.2. 
The frustration is no longer negligible. As one can see, 
the use of an effective spin model to describe x(^) I0r 
a spin chain coupled to phonons is well justified. An 
easy way to fit experimental data is to use Pade approx- 
imants for high temperature series expansion instead of 
exact complete diagonalization. The maximum is well- 
described by this procedure, even though the limit of very 
low temperatures is problematic. 

Figures fj] and || show additional results for g = 1.0 J 
and g — 1.5J and for oj = 10J. As J/uj is small the 
flow equation approach works especially well and thus 
the agreement is even better compared with the uj = 4 J 
results. 

Three final remarks regarding the susceptibility results 
are in order. In Ref. the authors observed that with 
increasing spin-phonon coupling the susceptibility curves 
are shifted to higher temperatures. The main reason for 
this is the static term oc (bj + bi) in the original QMC 
Hamiltonian (|l5|). This is not visible in our way of de- 
picting x{T) using the parameterization (Q) as we shifted 
the phonons to make this term vanish. If this static effect 
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FIG. 7: Magnetic susceptibility for g = 1.0J and lo = 10J. 
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FIG. 8: Magnetic susceptibility for g — 1.5 J and us = 10 J. 

is accounted for by the phonon shift, the positions of the 
maxima of x(T) are nearly the same for all parameters 
g and u> presented here. Only if g/ui and thus the ef- 
fective frustration are very large a small shift due to the 
temperature dependent frustration becomes visible. 

The second remark concerns the applicability of the 
flow equation procedure presented here. As we have al- 
ready seen in Fig. |2[ the ratio J/ui must not become 
too large. Thus, we are only able to show susceptibilities 
with J/lu < 1/4. The main problem for a thermodynamic 
quantity is the coth(u//2T) term in the formula ( |l0| ) for 
the effective couplings because it diverges for T — > 00. To 
reach higher values of J/ui a more subtle flow equation 
approach would be necessary. But even the straightfor- 
ward flow equation results agree well with the QMC data 
for not too large J /to and g/u. 

The third remark concerns the temperature depen- 
dence of the effective couplings. This dependence is fairly 
small for the parameter values in Figs. which is due 



to the small values J/ to. Thus good fits are also obtained 
neglecting the temperature dependence J(T) = J in a 
model of static effectpe couplings. This is in agreement 
with previous resultstEl. But we expect that for larger 
values of J/lu the temperature dependence will be im- 
portant. Yet such parameters are beyond the scope of 
the present approach based on the leading orders. Hence 
further investigations of this issues are called for. 



V. SUMMARY 

We applied two methods to study the phase diagram 
and the magnetic susceptibility of the isotropic antiferro- 
magnetic spin 1/2 Heisenberg model coupled to Einstein 
phonons via a bond coupling: the flow equation method 
and a quantum Monte Carlo approach. 

The flow equation approach maps the full problem onto 
an effective magnetic problem, i. e. onto a frustrated spin 
chain. The calculation was performed in leading order 
in g/u and in the two leading orders in J/lu and is thus 
systematic to order g 2 /lu 2 and g 2 J/iv 3 . The expansion 
approach works especially well in the antiadiabatic limit. 

To apply the flow equation formalism, it is important 
to use the spin operator which couples to the phonons 
in its normal-ordered form. To obtain x(T) from the 
calculated effective couplings Ji(T) and J2(T) we used 
high temperature series expansions and exact complete 
diagonalization. 

The quantum Monte Carlo calculations are completely 
independent of the effective couplings and are based on 
a loop algorithm. Spin and phonon degrees of freedom 
are dealt with an equal footing. 

The zero temperature phase diagram derived by these 
methods agrees well for intermediate J/lu and is in agree- 
ment with the result in Ref. |l4]. Approaching the limit 
J/lu —> we find a finite critical spin-phonon coupling 
g c /u as an important difference to the difference cou- 
pling. This is rigorous since a vanishing g c /tu for J/lu 
would correspond to the limit where our approach be- 
comes exact. We cannot fix the exact value of g c /LU in 
this limit as the flow equation approach corresponds to 
a leading order expansion in J/lu and g/u and thus is 
not exact for finite g/uu. Yet both our approaches yield 
a significantly lower value of g c /tu than the fourth order 
calculationlilJ. 

We compared the magnetic susceptibilities obtained by 
both methods and find good agreement. We conclude 
that at least in this parameter regime the thermodynamic 
properties of the Hamiltonian Eq. ^ can be described ex- 
cellently by an effective spin model with temperature de- 
pendent couplings. The restriction in g/u/ and J/ui for 
the the flow equations are due the leading order charac- 
ter of the approach presented here and could be overcome 
by a more advanced analysis. The flow equation method 
provides a good completion to an elaborate QMC analy- 
sis as it can be used more quickly, so that the fitting of 
parameters is easier. 
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